Irreversibility, Mechanical Entanglement and Thermal Melting in Superconducting 

Vortex Crystals with Point Impurities 
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We discuss the onset of irreversibility and entanglement of vortex lines in high T c superconductors 
due to point disorder and thermal fluctuations using a simplified cage model. A combination of Flory 
arguments, known results from directed polymers in random media, and a Lindemann criterion are 
used to estimate the field and temperature dependence of irreversibility, mechanical entanglement 
and thermal melting. The qualitative features of this dependence, including its nonmonotonicity 
when disorder is sufficiently strong, are in good agreement with recent experiments. 



Introduction 



The physics of vortex matter in high temperature su- 
perconductors (HTSC) has generated considerable inter- 
est since these materials were discovered about a decade 
ago jjj. The magnetic flux lines in a clean material 
prefer to sit on the sites of a triangular lattice forming 
an Abrikosov solid below the superconducting transition 
temperature T c . However, the combined effects of ther- 
mal fluctuations and various types of disorder, such as 
oxygen vacancies, columnar pins or twin boundaries, sig- 
nificantly alters the phase diagram, which is further com- 
plicated by typically large crystallographic anisotropics 
In fact, pinning of flux lines by disorder is important 
technologically, since the flux lines must be immobilized 
in order to eliminate dissipative losses associated with 
flux line motion. In this paper, we will focus on pinning 
by pointlike defects, such as oxygen vacancies. There are 
many recent studies in the literature that have considered 
other types of disorder HQ. 

It is well known that even an arbitrarily small amount 
of point disorder is sufficient to disrupt the delicate long 
range translational order of the Abrikosov crystal ||. 
However, the existence of a distinct "vortex glass" phase 
HQ at low temperatures (one which presumably con- 
tains dislocations and/or disclination defects) is still con- 
troversial. Indeed, experiments in which electron irradi- 
ation injects point disorder in twin-free YBa2Cu307_a 
(YBCO) samples in amounts large enough to destroy the 
first order melting transition have failed to find a sharp 
vortex glass transition with universal exponents || . On 
the other hand, the presence of a stable, topologically or- 
dered, i.e., dislocation-free, "Bragg glass" phase was re- 
cently suggested by various authors ||] at low magnetic 
fields |rp[ and temperatures in the presence of weak point 
disorder. We will loosely refer to such a configuration as 
the "ordered phase" . When dislocations are absent, the 
interaction of a flux line (FL) with its neighbors can be 
approximated by an effective confinement potential that 
keeps the FL from wandering away, similar to the Ein- 
stein model of phonons in conventional crystals. We ap- 



ply the Lindemann criterion to this "cage model" (where 
the confinement potential is approximated by a simple 
harmonic well) in order to investigate the combined ef- 
fect of thermal fluctuations and point disorder. A simi- 
lar model was used previously to investigate the effects 
of correlated disorder due to columnar pins in Ref. [[|. 
The discussion will focus primarily on Bi2Sr2CaCu20s+5 
(BSCCO) crystals, where the large anisotropy confines 
transitions out of the hypothetical Bragg glass to easily 
accessible fields and temperatures. We seek to under- 
stand the following three phenomena observed by exper- 
iments (See Fig. |l|): 
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FIG. 1. Schematic phase diagram indicating the approxi- 
mate location and shape of the irreversibility line (dashed) 
and the combined melting-entanglement line (solid) . The ar- 
row indicates the path taken for the situation depicted in 
Fig. | 

(i) The melting line (ML) marks the loss of order in 
the FL system, and is widely believed to be a first-order 
thermodynamic phase transition, most notably marked 
by a finite jump in the magnetization as a function of 
both temperature and applied field pi| . 

(ii) The irreversibility line (IL) marks the cross-over 
when a sample falls out of equilibrium on laboratory time 
scales. It is typically determined by the onset of hystere- 
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sis in the magnetization curve, and is recently demon- 
strated to be separate from the melting line, provided 
that one excludes effects due to surface barriers [ p"2"| |. 

(Hi) The entanglement line (EL) has only been stud- 
ied more recently and can be defined experimentally by 
a sudden increase in the critical current when the field is 
increased at low temperature JT^ | . Fairly independent of 
temperature, its driving mechanism may be the gradual 
proliferation of large dislocation loops favored by point 
disorder [j9|]l4fl . The onset of the large critical currents at 
low temperatures with increasing field coincides approxi- 
mately with the field above which the intensity of Bragg 
peaks in neutron diffraction experiments suddenly de- 
crease [Tq] . Unlike first order melting, these phenomena 
are not accompanied by a magnetization jump, suggest- 
ing a continuous phase transition or a crossover |b| . The 
proliferation of large scale dislocation loops would gener- 
ate entanglement via their screw dislocation components 
(see Fig . ||), as suggested years ago by Wordenweber and 
Kes 16 1 and by Brandt M| in the context of dimensional 
crossover. Provided line crossing barriers are large, such 
long wavelength entanglements would increase the effec- 
tiveness of pinning, similar to ideas about enhanced pin- 
ning in arrays of splayed columnar defects pl|. As the 
dislocations become denser with further increases in the 
field, vortex crossing angles will become large, flux cut- 
ting easier and critical currents should drop. The critical 
current in experiments does in fact gradually decrease 
with further increases in magnetic field [p"3| , indicative of 
a gradual evolution to a decoupled regime of "supercntan- 
gled" pancake vortices. 




Top view 
of cross section 



Side view 



FIG. 2. Dislocation loop in a topologically ordered vor- 
tex crystal. Nucleation of dislocation pairs in a constant z 
cross-section (corresponding to the bottom of the loop) be- 
comes more likely when the Lindemann criterion is reached. 
Dashed lines represent flux lines just behind the plane of the 
figure. Top view shows the Burgers' construction applied to 
a cross-section through the loop, dimensions. 

Numerical evidence for this mechanical entanglement 
scenario has been provided by simulations of Gingras and 
Huse of a disordered XY model jl4j, where vortex rings 
play the role of dislocation loops in an Abrikosov lattice. 
These authors find evidence for a continuous phase tran- 
sition to a defect free phase at low temperatures below 



a finite disorder strength in three dimensions. In d = 2, 
the transition occurs only at vanishing disorder strength. 
To make connections with real vortex arrays, effects of 
disorder in a given sample are assumed to increase with 
increasing magnetic fields. Direct simulations of vortex 
arrays with point disorder in three dimensions by Ryu et 
al. p~9| also support this picture, finding evidence for a 
proliferation of lines of disclination pairs (i. e., disloca- 
tions) above a critical filed threshold j2(| • 

Our model suggests a relatively simple scenario in 
which the melting, irreversibility and entanglement lines 
can be understood and connected in a way consistent 
with experimental observations. In this scenario, a sud- 
den increase in intravalley energy barriers is responsible 
for the portion of the IL that lies within the ordered 
phase. The lines ML and EL in Fig. |l| are connected, 
and both signal the proliferation of dislocations. How- 
ever, the primary mechanisms are different. The ML oc- 
curs at temperatures above the IL where point disorder 
is ineffective at pinning FLs. Melting is thus a competi- 
tion between thermal fluctuations and interaction effects. 
The thermal melting line is only slightly perturbed from 
its position in the absence of disorder. However, at tem- 
peratures below the IL, disorder-induced wandering of 
the FL beyond the Lindemann criterion possibly causes 
the proliferation of frozen-in dislocations that give rise 
to a mechanically entangled configuration with a sub- 
stantially increased critical current. This phenomenon, 
whether a crossover or a sharp continuous phase transi- 
tion, results from a competition between interactions and 
the effects of point disorder. 

For simplicity let us assume that the average magnetic 
field B is along the c-axis of the sample. The Hamilto- 
nian of a single FL whose configuration r(z) fluctuates 
about the line r(z) — rg can be written approximately as 



H[t{z)] = J dz 



dr 

2 \ dz 



+ V(r(z),z) 



r-r 



(1) 



where e = eo/7 2 is the line tension of the FL, — 
m z /m± >• 1 is the mass anisotropy, e = (^o/^ttA) 2 , 
and fc « eo/ag is the effective spring constant due to in- 
teractions with neighboring FLs, which are typically at a 
distance eto = \J <&q/B. This approximation to k requires 
B > H c i , so that the vortices interact with a logarithmic 
potential at the nearest neighbor spacing. The proba- 
bility of a particular vortex trajectory is proportinal to 
exp{— 7i[r(z)]/fcsT}. (The Boltzmann constant ks will 
be suppressed henceforth.) The disorder potential V is 
assumed to have zero mean and short-range correlations 



(V(r,z)V(r',z')) = ASf \r - r')S{z - z'), 



(2) 



where 5^ (r) is a delta function smeared out to the vortex 
core diameter £. It will be useful to define a dimensionless 
disorder parameter given by 
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A = A/(e^). 



(3) 



In the absence of the cage potential (k = 0), the model 
reduces to a directed polymer in a 2+1 dimensional ran- 
dom medium, which has been studied extensively |2l) . At 
low temperatures, the FL makes transverse excursions 

u 2 {£) ee (\r(z) r(z + £)\ 2 ) a £ 2 (44) 2C , I > 4, (4) 

and the energy barriers between low-energy configura- 
tions grow as |22[ 



Upo(£)f*T*(t/t, 



\2C-1 



(5) 



A Flory-type argument |23| , combined with the known 
value of C ~ 5/8 pi) gives 



4 » (£/ 7 )(A 7 )- 1/3 , 
T* a (^o/7)(A 7 ) 1/3 . 



(6) 
(7) 



At shorter length scales, where it is justified to expand 
the random potential in Eq.(|l|) about r = ro, 

V(v(z), z) « V(r , z) + (r(z) - r ) • V r V(r = r , z), (8) 

the fluctuations are given by the Larkin-Ovchinnikov for- 
mula 



(9) 



However, this regime is not relevant to Lindemann dis- 
cussions of entanglement provided that c^ao > £, and 
will be ignored henceforth, (cz, is the Lindemann con- 
stant, see below.) Thermal noise is ineffective at de- 
pinning the FL from individual pinning centers as long 
as the magnitude of thermal fluctuations are smaller 
than quenched fluctuations at the crossover length scale, 
i.e. when u 2 h (£ c ) — T£ c /i <C £ 2 , or equivalently when 
T <C T* . At higher temperatures, the excursions of the 
FL at short distances are determined mainly by thermal 
noise. However, the scaling of transverse fluctuations is 
still determined by point disorder when £ > £ C (T), where 
the temperature dependent crossover length scale is [9 



4(T) = 



£ 11 e c(T/T>) 



T < T*, 
T* < T, 



(10) 



where c is a constant of order unity. The sharp exponen- 
tial increase of £ C (T) as T exceeds T* is caused by ther- 
mal smearing of the disorder potential at the marginal 
dimension d = 3. Transverse fluctuations of the FL are 



](£) 



T >T* 



£>lc 

i < 4(T), 



tc(T) 
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T > T*, £ > £ C (T), 



(11) 



and typical energy barriers are 



U p0 (£) 




T < T*, I > 4, 

T > T*, i < 4(T), (12) 

T > T*, I > 4(T). 



The main effect of the cage potential is to block trans- 
verse excursions of the FL beyond a confinement length 
£* a \Jejk = ao/7, which is obtained by balancing the 
tilting and confinement energies. Every time the vortex 
wanders a distance I* along the z— axis, it is reflected 
back by the walls of the cage potential and must restart 
its thermal or disorder dominated random walk. Thus, 
the mean square displacement of the FL around its equi- 
librium position and the typical energy barriers are re- 
spectively given by 



U P (T) 



:ul{t), 

■u p0 {t) 



(13) 
(14) 



The vortex line acts as if it were broken up into inde- 
pendent subsystems of length £* . Similar arguments for 
the effect of interactions on point disorder and thermal 
fluctuations were given for flux liquids in Ref. [p4[. 



The Irreversibility Line 

Using the results ( ^OfjT^ ), let us investigate the phe- 
nomena described earlier in the context of this simple 
model. Naively, one would expect the IL to be associ- 
ated with a sudden increase in energy barriers against 
the motion of individual FLs. When T where 
T^ 1 ) is defined by the condition 



4(T« 



£* 



(15) 



pinning centers do not offer a gain in the free energy of 
the FLs, and there are no barriers against FL motion. 
The slowest relaxation time of the FL in this regime is 



?7 07TA ag 



(16) 



where 7/ = $0^2 /(PnC 2 ) [£5[ is the friction coefficient 
of the FL, and p n is the resistivity of the normal elec- 
trons in vortex cores. When T drops into the range 
T* < T < T^ 1 ' , energy barriers grow rapidly and the 
relaxation time increases to 



r(T) = r e' 



= ^M P (T)/T 



tq exp ■ 



7 c (t( 1 ); 



1/4' 



(17) 

T*<T<T {1 \ (18) 



Eventually, energy barriers saturate at T sa T* , resulting 
in 



T* ft 

r(T)=roexp - - 



1/4- 



(19) 
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Thus, in a range of temperatures near T*, the relaxation 
time increases rapidly to macroscopically long times. The 
irreversibility condition is satisfied when log(r(T)/ro) 
becomes fairly large. Upon further manipulation, this 
yields 



Bil{T) 



B IL {T*)exp{-2c{T/T*) 3 }, 
B IL {T*)(T/T*)' 8 , 



T > T*, 

rp ^ rp* 



(20) 



The Melting and Mechanical Entanglement Lines 

In order to estimate where dislocations appear in the 
B — T phase diagram, we can employ a Lindemann cri- 
terion as follows: 



u 2 (T m (B))=c 



2 2 
L"0' 



(21) 



where cl ~ 0.15 — 0.2 is the phenomenological Linde- 
mann constant J2(|. Thus, for a given sample, the form 
of u 2 (T) will determine the shape of the melting line. 
We next discuss various regimes in this function and the 
corresponding melting fields B m {T). 

1. T <C T*: In this region, the Lindemann criterion 
(pl|) becomes 



Upon simple manipulation, this yields 



10/3 



„ 16 / 3 



(22) 



(23) 



Taking, as parameters for BSCCO, f = 20 A, 7 = 100, 
e £ = 1000 K, and T* = 10 K, c L « 0.17 leads to 
S m (0) « 400 Gauss, in reasonable agreement with ex- 
perimental results jl3j . Note that the results depend sen- 
sitively on the value of the roughness exponent ( k 5/8 
and the Lindemann constant cl- 

2. T* < T < TW: In this region, thermal fluctuations 
spread the vortex probability distribution and weaken the 
pinning effect of point disorder, but the relevant length 
scale at which the cage potential is seen is still determined 
by the anomalous wandering characteristics of point pin- 
ning. The result is 



B m (T) « B, 



10/3 



e ¥ (T/T*) 3 



T* <T < T (1) 



(24) 



indicating an increase in the melting field with T. 

3.T 3> T^- 1 ': In this region, the FL sees the cage po- 
tential before it has a chance to take advantage of point 
disorder, which is effectively smeared out and can be ne- 
glected as a first approximation. We then find the well- 
known portion of the melting line, 



B m (T) 



7 2 T 2 



(25) 



Note that eo has an implicit dependence on temperature, 
and the effective anisotropy parameter 7 can depend on 
the amount of point disorder. The melting field is a de- 
creasing function of temperature in this regime. 

Combining these three regimes gives rise to a non- 
monotonic melting field B m (T), which is qualitatively 
depicted in Fig. |l|. At the lower temperature side of 
this peak, B m (T) lies below the irreversibility line, and 
configurations of the FL are effectively pinned by point 
disorder. Thus, it is more appropriate to regard the Lin- 
demann criterion in this regime as a pointer to the on- 
set of mechanical entanglement, rather than melting p7[ . 
More sophisticated methods are needed to determine the 
nature (first or second order transition, or crossover) of 
this onset. The high temperature side of the curve corre- 
sponds to the more familiar and studied phenomenon of 
flux lattice melting, which was observed as a first-order 
thermodynamic phase transition. Finally, the nonmono- 
tonic behavior disappears when the disorder is too weak, 
i.e. for T* < cl^co/7- The transition line observed re- 
cently by Khaykovich et al. in highly anisotropic Bi- 
bascd HTSCs is in good qualitative agreement with our 
model. 



Numerical Results 



In order to test the analytical results, we have numer- 
ically determined u 2 (T) for various values of the mag- 
netic field. The general procedure is outlined below, a 
more detailed discussion on Eqs. ( p6| - p^ ) can be found in 
Ref. [§§. 

The "weight" function for a FL starting from at 
z = and wandering across a sample of thickness £ to 
position r is given by 

r r'(t)=r 

W(r,Ti;£)= Vr'(z)exp{-n[r'(z)]/T}, (26) 

Jr'(0)= ri 

which satisfies a Schrodinger-like equation with a time- 
dependent ("time = z) potential, 



d z W{r,Yi;z) = 



T ^2 k , 

— V r 

2e r 2 1 



V(t(z),z) 



W^Vi-z). 
(27) 



For a given realization of point disorder, the weight func- 
tion can thus be determined by discretizing and step- 
ping this equation forward in the time-like parameter z. 
The probability distribution for the position of a FL at 
a height z which is far from the boundaries of a sample 
of thickness £ is proportional to 



4 



P(r; z) oc / dr, J dv f W{v f , r; £ - «)W(r, r 4 ; z). (28) 

For a given realization of disorder, the weight function 
ratio W(r,Ti,z)/W(r' ,Ti; z) becomes independent of 
when z>f since the FL forgets its initial position after 
bouncing off of the confinement potential a few times. 
We have numerically confirmed that the starting posi- 
tion is indeed indistinguishable for z > 5£*, in agreement 
with the argument above. Thus, the integral over ini- 
tial and final positions in Eq. ( pa ) can be eliminated by 
setting Yi = Yf = Yq. (Alternatively, one could start 
from a uniform weight function.) Furthermore, the two 
weight functions that appear in Eq. (|2|) are statistically 
independent, since they traverse different portions of the 
sample. Taking advantage of these properties, we have 
determined P by multiplying two weight functions which 
were propagated through different samples by z = 5£*. 
This probability distribution was used to determine the 
fluctuation amplitude, which was subsequently averaged 
over different realizations of disorder. Fig. || shows the 
results for the parameter values A — 7 — 1, for different 
values of the dimensionless field b = B^ 2 /$o- In order to 
solve Eq. (p^), each z— slice was discretized to a square 
lattice with lattice constant £, and consecutive z— slices 
were a distance Sz — 0.05 £/7 apart. The random poten- 
tial V was determined independently at each site of the 
square lattice and had a correlation length of £/7 along 
the z— direction. 

0.30 1 . 1 -7 1 . n 



over many more samples to achieve reliable results. Nev- 
ertheless, reentrant behavior can be deduced from the 
clear nonmonotonicity of the fluctuation amplitude as a 
function of temperature, as seen in Fig. |3[ This feature 
can be understood qualitatively as follows (See Fig. ||): 
As the temperature is increased from T = 0, thermal 
fluctuations will initially act to blur the pinning centers, 
thus reducing their pinning strength and giving rise to a 
straighter FL with a smaller fluctuation amplitude. How- 
ever, further increases in temperature eventually spread 
the FL, increasing the fluctuation amplitude. 
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FIG. 4. Temperature scans of the 10 samples (thin lines) 
used to determine the disorder averaged fluctuation ampli- 
tude (thick line) for b — 1/900. The spread in data points 
increases significantly when T drops below T*. 
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FIG. 3. Results of a transfer matrix computation of trans- 
verse fluctuations in the cage model. The curves represent 
temperature scans at b = _B£ 2 /<E>o = 1/100 (circles), 1/400 
(squares), 1/900 (diamonds), 1/1600 (upwards triangles), and 
1/2500 (downwards triangles). The parameters in this run 
were A = 7 = 1. cl = 0.15 is shown as a representative 
Lindemann constant. 

While this method allows exact thermal averaging of 
fluctuations, obtaining the disorder average is much more 
difficult, as illustrated in Fig. ^. The spread of data 
points (each representing one realization of disorder) in- 
creases dramatically for T < T*, when disorder becomes 
dominant. Ideally, disorder averaging should be done 
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FIG. 5. Transverse fluctuations of a FL are minimized at 
an intermediate temperature (middle) due to the competition 
between point disorder and thermal noise. Reentrant behav- 
ior is possible if the Lindemann criterion is crossed during a 
temperature sweep (See Fig. ^). 



Discussion 

In this paper, we have presented the results of a sim- 
ple "cage model" for a FL embedded in a dislocation- 
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free vortex phase. The model describes the rela- 
tive motion of a vortex with respect to the cage 
provided by its neighbors, and is thus insensitive to 
the long wavelength disorder fluctuations studied in 
Refs. ■ An equivalent theory could presumably 

be constructed within the usual Debyc phonon model by 
setting (|u(r + a ei) — u(r)| 2 } = c|<Zq, where ei is a lat- 
tice vector and the brackets and overbar represent ther- 
mal and disorder averages respectively Q . We have dis- 
cussed the sharp increase in energy barriers (suggestive 
of an irreversibility line), and the onset of instabilities in 
the ordered phase according to a Lindemann criterion. 
One shortcoming of this model is its inability to verify 
the stability of the ordered phase with respect to large- 
wavelength fluctuations caused by point disorder. How- 
ever, the stability of the Bragg Glass phase with respect 
to such fluctuations has recently been demonstrated 
Thus, this model should predict approximately when a 
given sample falls out of equilibrium, and when disloca- 
tions proliferate, although more sophisticated models are 
necessary to determine whether the dislocation prolifer- 
ation occurs through a phase transition or a crossover. 

The low temperature entanglement field B m (0) should 
not be confused with the so called "decoupling field" 
H x w $ /(d 2 7 2 ) for highly anisotropic HTSCs, above 
which a discretized version of Eqs. ( pq ) and (27) should 
be used |2S| ]. vortex line behave roughly independently 
in each Cu-0 layer. The ratio of these two fields is 
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For materials such as YBCO or BSCCO, this ratio is 
typically less than 1. A more detailed account of layer- 
ing effects and the dispersive nature of the line tension is 
beyond the scope of this work. 
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